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ABSTRACT 

We discuss how inhomogeneities of the universe affect observations of the 
gravitational lensing; (1) the bending angle, (2) the lensing statistics and 
(3) the time delay. In order to take account of the inhomogeneities, the 
Dyer-Roeder distance is used, which includes a parameter representing the 
dumpiness of the matter along the line of sight. It is shown analytically that 
all three combinations of distances appearing in the above observations (l)-(3) 
are monotonic with respect to the dumpiness in general, for any given set of 
the density parameter, cosmological constant and redshifts of the lens and 
the source. Some implications of this result for the observation are presented; 
the dumpiness decreases both the bending angle and the lensing event rate, 
while it increases the time delay. We also discuss cosmological tests using the 
gravitational lensing in the clumpy universe. 
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Introduction 



It is one of the most important longstanding problems to determine the cosmological 
parameters (Weinberg 1972; Peebles 1993). There are some methods to determine the 
cosmological parameters by using the gravitational lenses (Refsdal 1964a, 1964b, 1966; 
Press and Gunn 1973; Blandford and Narayan 1986; Fukugita, Futamase and Kasai 1990; 
Fukugita, Futamase, Kasai and Turner 1992; Schneider, Ehlers and Falco 1992; hereafter 
SEF). Most of them can be classified into the following three typical observations: (1) the 
bending angle, (2) the lensing statistics and (3) the time delay. It is of great importance to 
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clarify the relation between the observation in the realistic universe and the determination 
of the cosmological parameters. In particular, it has been discussed that inhomogeneities 
of the universe may affect the cosmological tests (Dyer and Roeder 1972, 1974; hereafter 
DR; Schneider and Weiss 1988a, 1988b; Linder 1988; Futamase and Sasaki 1989; Kasai, 
Futamase and Takahara 1990; Bartelmann and Schneider 1991; Watanabe, Sasaki and 
Tomita 1992; Sasaki 1993). In this paper, we use the Dyer- Roeder (DR) angular diameter 
distance in order to take account of the inhomogeneities (DR). Besides the three parameters 
(the Hubble constant, the density parameter and the cosmological constant), the DR 
distance contains an extra parameter representing the dumpiness of the matter along the 
line of sight. Two extreme cases can be represented by the DR distance; one is the distance 
in the Friedmann-Lemaitre-Robertson- Walker (FLRW) universe, the so-called filled beam, 
while the other is that of the so-called empty beam when the right ray propagates through 
the empty region. For comparison with the filled beam, the empty beam has been frequently 
and numerically used in the literature (For instance, Fukugita, Futamase and Kasai 1990; 
Fukugita, Futamase, Kasai and Turner 1992). However, it has not been clarified whether 
the observed quantities and/or the cosmological parameters obtained in the arbitrary case 
of the dumpiness parameter are bounded between those for the filled beam and the empty 
beam. Moreover, some cases of dumpiness parameter have been investigated numerically 
for fixed redshifts of the lens and the source (For instance, Alcock and Anderson 1985). 
However, the effect of the dumpiness on the observable depends on the redshifts of the lens 
and the source. Therefore, it is important to clarify how the observation of gravitational 
lensing depends on all the parameters (the density parameter, cosmological constant, 
dumpiness parameter and redshifts of the lens and the source), since the dependence 
seems complicated. Hence, we investigate analytically the arbitrary case of the dumpiness 
parameter, for any set of the density parameter, cosmological constant and redshifts of the 
lens and the source. 

This paper is organized as follows. Section 2 introduces three types of distance 
combinations appearing in the gravitational lensing (l)-(3). Section 3 shows the basic 
equations for the DR distance. Section 4 presents the proof that all these distance 
combinations are monotonic functions of the dumpiness parameter. We also clarify the 
effect of the dumpiness on the observations (l)-(3). In section 5, we discuss how the 
dumpiness affects the determination of the cosmological parameters. Conclusions are 
summarized in section 6. 
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2. Distance combinations appearing in gravitational lenses 

(1) bending angle 

The lens equation is written as (SEF) 



Here, (3 and 6 are the angular position vectors of the source and image, respectively, and 
a is the vector representing the deflection angle. Dls and Dos are the angular diameter 
distances from the lens to the source, and from the observer to the source, respectively. The 
effective bending angle (D LS / Dos)ot appears when we discuss the observation concerning 
the angle such as the image separation and the location of the critical line (Blandford and 
Narayan 1986, SEF). Hence the ratio Dls I 'Dos plays an important role in the discussion on 
the observation concerning the angle. It has been argued that, in calculating the bending 
angle, the density along the line of sight should be subtracted from the density of the lens 
object (Sasaki 1993). However, we assume that the density of the lens is much larger than 
that along the line of sight, so that the effect of the dumpiness on a can be ignored. Thus, 
we consider only the ratio Dls /Dos in the following. 

(2) lensing statistics 

The differential probability of lensing events is (Press and Gunn 1973; SEF) 



where ul is the number density of the lens, dl is the physical length of the depth and a is 
the cross section proportional to DolDlsI Dos- Here, Dol is the angular diameter distance 
from the observer to the lens. Since dl depends only on the cosmological parameters in the 
FLRW universe, we investigate the combination D Q lDls/ Dos m order to take account of 
the dumpiness of the matter. 

(3) time delay 

The time delay between two images A and B is written as (Refsdal 1964b, SEF) 



dr = ariLdl, 



(2) 



&t AB = 



1 + z L DqlDqs 
c D LS 




— ct 



(3) 



where at a and olb are the bending angles at the images A and B, respectively. 
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3. DR distances 

The DR angular diameter distance is determined by (DR, SEF, Sasaki 1992) 

d 2 - 3. . 5 



(Ji(ft D + -{l + z)>a&W = 0, (4) 

where the parameter a represents the dumpiness of the matter along the light ray. In the 
FLRW universe, a is unity, while a vanishes when the light ray propagates through the 
empty space. Here w is an affine parameter, which is assumed to be that in the FLRW 
universe, determined by 

^ = (1 + z) 2 y/nz(l + zf - \z{2 + z) + (1 + z)\ (5) 

where Q and A denote the density parameter and the cosmological constant, respectively. 
Since the coefficient of the last term of Eq.([|), 3aQ/2, comes from the Ricci focusing by the 
matter along the line of sight, the DR angular diameter distance is a decreasing function of 
a for a fixed redshift (DR, SEF). That is to say, 

-Dol(«i) > A3l(«2) for ai < a 2 . (6) 

In reality, the parameter a takes various values according to mass distribution. For 
instance, it can take a rather low value such as 0.5 in the clump model (Kasai, Futamase 
and Takahara, 1990; Linder, 1997). 

The DR equation @ must be solved under the boundary conditions, 

D(z 1 ,z l ) = 0, (7) 



and 



—D(zi, z 2 ] 
dz 2 



a(zi) 



H{z x ) 



where a(zi) and H(z\) denote the scale factor and the Hubble constant at z\, respectively. 



4. Monotonic properties 

(1) D LS /D os 

It has been shown that the distance ratio Dls/Dos satisfies (Asada 1997) 

^^-{otx) < ^^-{a 2 ) for ai < a 2 . (9) 
^os ^os 
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This is proved as follows. For fixed zs, £1 and A, the ratio D LS /D s can be considered as a 
function of z L , X a (zi). We define Y a {zi) as Dsl/Dos, where D S l is the DR distance from 
the source to the lens. Owing to the reciprocity (SEF), we obtain 

Y a (z L ) = l±?*X a (z L ). (10) 

Since Y a depends on Zl only through D$l, it obeys the DR equation, 

<r Y a (z L ) + 1(1 + z L ) 5 aQY a (z L ) = 0, (11) 



where w l is an affine parameter at the lens. We define the Wronskian as 

^..nj = (y M f*-r,§*). (12) 



Using Eq . ([TTD , we obtain 

d 



■W(Y ai ,Y a2 ) < for ai <a 2 . (13) 



dw L 

Since both Y ai and Y a2 vanish at Zl = zs, we obtain 

W(Y ai (z s ),Y a2 (z s )) = 0. (14) 

From Eqs.(|T3"D and (H), we find 

w(y ai ,y aa ) >o, (is) 

where we used the fact that the affine parameter w defined by Eq.fl5|) is an increasing 
function of z. Equation ( |15D is rewritten as 

(/ ln^>0. (16) 



dw L Y ai 

Since Y a always becomes 1 + zs at the observer, we find 



From Eqs.flTED and (|TTD, we obtain 

In ^ > 0. 



Y 



This leads to 

1 ' > I- (19) 



L «2 
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where we used Eq.([H]). Thus, Eq.(|9]) is proved. 

It should be noted that Eq. flT8"| ) holds even if one uses the opposite sign in the definition 
of the affine parameter in Eq.(||). From Eqs.(|lD and (||), the image separation as well as the 
effective bending angle increases with a. 

(2) DolD L s/Dos 

Next let us prove that DolDls I Dos increases monotonically with a. We fix Q, A, zl 
and zs- Then it is crucial to notice that the distance from the lens to the source can be 
expressed in terms of the distance function from the observer, D(z), as (Linder 1988) 

c r Ws dw 

Dls = tt(1 + z L )D 0L D s / (20) 
H Jw L D(z) 2 

where Hq is the Hubble constant at present. It is verified in a straightforward manner that 
Eq.(^0|) satisfies the DR equation (f|) with the boundary conditions (0) and (||). Equation 
( pOl) is rewritten as 

Dos -"0 Jw L V(z) 2 

The right hand side of this equation depends on a only through Dol/ D(z). By the similar 
manner to the proof of Eq.(|9|), we obtain for Zl < z < z$ 

D{z) D(z) 
By applying Eq.(p2|) to Eq.(|2l|), we obtain 

DolDls , , D Q lD L s , , f ^ / o0 x 

^OS u OS 

Therefore, the gravitational lensing event rate increases with a. 

(3) Dls/DolDos 

Finally, we investigate the combination of distances appearing in the time delay. By 
dividing Eq.(fD by Eq. @, we obtain 

DolDos , \ ^ DolDos , \ f ^ , nA \ 

— Fi ("i) > — Fi ("2) for ai < a 2 . (24) 

L>LS L>LS 

Thus, the time delay decreases with a. 

Before closing this section, a remark is given: As shown above, the three types 
of combinations of distances are monotonic with respect to the dumpiness parameter. 
However, some of other combinations of distances are not monotonic functions of a, though 
these combinations could not be necessarily relevant to physical problems. For instance, a 



combination D L sl y cDqs/ H is not a monotonic function of a. 
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5. Implications for cosmological tests 

In this section, we consider the three types of the cosmological test which use 
combinations of distances appearing in the gravitational lensing. Let us fix the density 
parameter in order to discuss constraints on the cosmological constant. 

(1) D L s/Do S 

The following relation holds 

^(Ax)<^(A 2 ) for Ai < A 2 . (25) 
This is proved as follows. Let us define 

D LS (a,fi } A) 

X ^ z l) = n — 7 — pTTT' ( 26 ) 
Dos{oi,il, A) 



and 



D fir(a,S2, A) 



By the reciprocity (SEF), we obtain 



W = ^^a(^). (28) 

J- T 



The ratio Ia^l) satisfies 

r/ J 3 



F A (^ L ) + -(1 + z L ) b aQY a (z L ) = 0. (29) 



dw 2 L ~" y ' 2 



For Aj (i = 1, 2), the affine parameter Wi satisfies 



dw, 
We obtain 



(llL ~ (1 + z L f Jnz L (l + z L y - X lZL {2 + z L ) + (1 + z L f. (30) 



du>l du> 2 

We define the Wronskian as 



d ^ > p. f 0r A, < A 2 . (31) 



Using Eq . ([29]) and Eq . (j3ll) , we obtain 

d 



dzi 



W(Y Xl ,Y X2 )<0 for A!<A 2 . (33) 
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Since Y x always vanishes at Zl = z$, we obtain 

W(Y Xl (z s ),Y X2 (z s )) = 0. 

From Eqs.(P5|) and (P4|), we find 

H/(y Al) y A J>0 for A!<A 2 . 



This is rewritten as 



where we used Eq.( 



d ln^>0 for A!<A 2 , 



dz L Y Xl 

Since Y x always becomes 1 + z$ at the observer, we find 



ln*f' = °>=0. 

Y Xl (z L = 0) 



(34) 

(35) 
(36) 

(37) 



From Eqs.fl36"D and fl37|), we obtain 



Y) 



In ^ > for Ai < A 2 . 



This leads to 



X 



A 2 



> 1 for Ai < A 2 



where we used Eq.(|^). Thus, Eq.fl25|) is proved. 



(38) 



(39) 



Together with Eq.(|^), Eq. (|25|) means that, in the cosmological test using the bending 
angle, the cosmological constant estimated by the use of the distance formula in the FLRW 
universe is always less than that by the use of the DR distance (0 < a < 1). 

(2) DolD L s/Dos 

By multiplying Eq.(p5|) with 

Dol(Xi) < A?l(A 2 ) for A!<A 2 , (40) 

we obtain 



DolDls / x x ^ D OL D LS 

-(Ai) < — (A 2 ) for Ai < A 2 . 



Dos D os 

Equation ( f40| ) can be proved, for instance in the following manner: The DR distance is 
written as the integral equation (Schneider and Weiss 1988a, Linder 1988) 



(41) 



3 c f z dw 

D(z;a) = D(z;a=l) + ~—(l-a)n dy — 

2 Hn Jo dz 



(l + y) i D(y,z-a = l)D(y 1 a). (42) 



z=y 



- 9- 



This is rewritten as (Schneider and Weiss 1988a, Linder 1988 

00 T3 c 
D(z; a) = D(z; a = 1) + ^ 

i=l 

where Ki(y,z) is defined as 



2 K 



o 



;i - a )n 



f dyKi(y,z)D(y;a = 1) 

JO 



(43) 



;i + x) 4 J D(x,y;a = 1) 



and 



K i+ i(x,y) = dzK 1 (x,z)K i (z,y). 

J X 

From Eqs.fl3~iP, (|44]) and ([4^), it is shown that for x < ?/ 

Ki(x,y; Ai) < Ki(x,y; A 2 ) for A x < A 2 , 
where we used the following relation in the FLRW universe as 



(44) 



(45) 



(46) 



D(x, y;a—l, Ai) < D(x, y;a — l, A 2 ) for Ai < A 2 . (47) 

Using Eqs.(H), (0) and (0), and the positivity of i^, we find Eq.(gg). 

From Eqs.(p3|) and (^H), it is found that, in the cosmological test using the lensing 
events rate, the cosmological constant is always underestimated by the use of the distance 
formula in the FLRW universe. 



(3) D L s/DolD 



OS 



We consider the combination D LS / D OL D Q s- Since the time delay is measured and the 
lens object is observed, DqlDosI 'D ls can be determined from Eq.(|3|). On the other hand, 
when we denote the dimensionless distance between z\ and z 2 as d\2 = HqD^/c, which does 
not depend on the Hubble constant, we obtain 



DqlDqs c dohdos 



D 



LS 



Hq d LS 



Then, Eq.(^) becomes 



dohdos , x ^ dohdos i x c 

-(cti) > — : (a 2 ) for a\ < a 2 . 



d 



LS 



d 



(48) 



(49) 



LS 



Thus, from Eqs.(|48|) and fl49|), it is found that H estimated by using the DR distance 
decreases with a. Thus, the Hubble constant can be bounded from below when we have 
little knowledge on the dumpiness of the universe. The lower bound is given by the use of 
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the distance in the FLRW universe. On the other hand, the combination D LS /D 0L D s is 
not a monotonic function of the cosmo logical constant. Therefore, the relation between the 
dumpiness of the universe and the cosmological constant is not simple, since it depends on 
many parameters (zl, z$, fi, A, a). 

It should be noted that even the assumption of the spatially flat universe (Q + A = 1) 
does not change the above implications for the three types of cosmological tests, since the 
cosmological constant affects the DR distance formula only through the relation between z 
and w, Eq.fl^). 



6. Conclusion 

We have investigated the effect of the dumpiness of the matter along the line of sight 
on the three types of observations on (1) the bending angle, (2) the lensing statistics 
and (3) the time delay. First, it has been shown analytically that the combinations of 
distances which appear in the gravitational lensing (l)-(3) are monotonic with respect to 
the dumpiness parameter a when the DR distance is used. This property presents us with 
the following implications for the observation: In the clumpy universe approximated by the 
DR distance, the bending angle is smaller, the lensing events occur less frequently and the 
time delay is longer than in the FLRW universe. In the cosmological tests using (1) the 
bending angle and (2) the lensing statistics, the use of the DR distance always leads to the 
overestimate of the cosmological constant. However, the same is not true of the time delay, 
since the combination of DR distances in the time delay is not monotonic with respect 
to the cosmological constant. Whether the effect of the dumpiness parameter enhance 
the effect of the cosmological parameter or not depends on other parameters. Rather, the 
primary cosmological use of the time delay is the physical estimation of the Hubble constant 
(Refsdal 1964b, Blandford and Kundic 1996). It has been found that the use of the DR 
distance never lowers the estimate of the Hubble constant from the time delay. 

We have taken the dumpiness parameter a as a constant along the line of sight. 
However, as an extension of the DR distance, a can be considered phenomenologically as 
a function of the redshift in order to take account of the growth of inhomogeneities of the 
universe (Linder 1988). However, in the above consideration, particularly in the proof of the 
monotonic properties, it has never been used that a is constant. Hence, all the monotonic 
properties and the implications for cosmological tests remain unchanged for variable a (z). 
That is to say, when a±(z) < a^iz) is always satisfied for < z < zg, all we must to do 
is to replace parameters ati and a<i with functions ol\(z) and ot^iz) in Eqs.@, (23) and 



In particular, when a(z) is always less than unity on the way from the source to the 
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observer, both of the three combinations of distances appearing in (1) and (2) are less, and 
the combination in (3) are larger than those in the FLRW universe. Then, the decrease in 
the bending angle and the lensing event rate, and the increase in the time delay hold even 
for a generalized DR distance with variable a(z). 

The DR distance is useful for theoretical and conceptual studies on the lensing in 
the clumpy universe. However, the DR distance seems too simple to describe the realistic 
universe where the density fluctuation is of stochastic nature. In fact, there are many lines of 
sight on which a may be deviating appreciably from unity (Kasai, Futamase and Takahara 
1990; Linder 1997). The conclusion obtained by the use of the DR distance encourages us 
to make a statistical discussion on the observation of gravitational lensing in the clumpy 
universe. For such a discussion, it is necessary to investigate the light propagation in the 
realistic universe numerically (Schneider and Weiss 1988a, 1988b; Bartelman and Schneider 
1991; Watanabe, Sasaki and Tomita 1992) and perturbatively (For instance, Seljak 1994; 
Bar-Kana 1996). 
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